Spectral analysis of finite-difference approximations of 1−d waves in non-uniform grids

Preserving the finite positive velocity of propagation of continuous solutions of wave equations is one of the key issues, when building numerical approximation schemes for control and inverse problems. And this is hard to achieve uniformly on all possible ranges of numerical solutions. In particular, high frequencies often generate spurious numerical solutions, behaving in a pathological manner and making the propagation properties of continuous solutions fail. The latter may lead to the divergence of the “most natural” approximation procedures for numerical control or identification problems. On the other hand, the velocity of propagation of high frequency numerical wavepackets, the so-called group velocity, is well known to be related to the spectral gap of the corresponding discrete spectra. Furthermore most numerical schemes in uniform meshes fail to preserve the uniform gap property and, consequently, do not share the propagation properties of continuous waves. However, recently, S. Ervedoza, A. Marica and the E. Zuazua have shown that, in 1− d, uniform propagation properties are ensured for finite-difference schemes in suitable non-uniform meshes behaving in a monotonic manner. The monotonicity of the mesh induces a preferred direction of propagation for the numerical waves. In this way, meshes that are suitably designed can ensure that all numerical waves reach the boundary in an uniform time, which is the key for the fulfillment of boundary controllability properties.